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DETERMINISTIC  NETWORK  INTERDICTION 

R,  Kevin  Wood 

Operations  Research  Department,  Naval  Postgraduate  School 
Monterey,  CA  93940,  U.S.A. 


Abstract — Interest  in  network  interdiction  has  been  rekindled  because  of  attempts  to  reduce  the 
flow  of  drugs  and  precmsor  chemicals  moving  through  river  and  road  networks  in  South  America. 

This  paper  considers  a  problem  in  which  an  enemy  attempts  to  maximize  flow  through  a  capacitated 
network  while  an  interdictor  tries  to  minimize  this  maximum  flow  by  interdicting  (stopping  flow 
on)  network  arcs  using  limited  resources.  This  problem  is  shown  to  be  NP-corapIete  even  wben 
the  interdiction  of  an  arc  requires  exactly  one  unit  of  resource.  New,  flexible,  integer  programming 
models  are  developed  for  the  problem  and  its  variations  and  valid  inequalities  and  a  reformulation  are 
derived  to  ti^ten  the  LP  relaxations  of  some  of  these  models.  A  small  computational  example  from 
the  literature  illustrates  a  hybrid  (partly  directed  and  partly  imdirected)  model  and  the  usefulness 
of  the  valid  inequalities  and  the  reformulation. 

1.  INTRODUCTION 

Interest  in  the  topic  of  network  interdiction  has  been  revived  recently  resulting  from  the  U.S. 
anti-drug  effort.  In  particular,  the  U.S.  Army’s  SOUTH  COM,  which  directs  U.S.  anti-drug  efforts 
in  South  America,  would  like  to  best  allocate  its  limited  resources  to  interdicting  coca,  partially 
processed  cocaine  and  precursor  chemicals  in  the  South  American  drug  producing  areas  [1]. 
Current  emphasis  is  on  interdicting  the  flow  of  precursor  chemicals.  Most  of  the  traffic  in  these 
chemicals  is  carried  along  rivers  in  remote  forested  regions  and  along  some  inter-connecting  roads. 
Consequently,  in  this  situation,  the  interdiction  problem  can  be  viewed  as  a  network  interdiction 
problem.  It  is  the  purpose  of  this  paper  to  review  earlier  deterministic  network  interdiction 
models,  devise  new  solution  techniques  for  these  models,  and  develop  new  models  and  solution 
techniques.  We  believe  that  these  models  may  be  useful  when  combined  with  simulations  for 
devising  randomized  strategies  for  interdicting  drug  and  precursor  chemical  traffic.  A  subsequent 
paper  will  investigate  probabilistic  and  game-theoretic  models. 

The  simplest  network  interdiction  problem  arises  from  the  well-known  max  flow-min  cut  theo¬ 
rem  [2].  In  this  case,  an  enemy  attempts  to  traverse  from  node  s  to  node  f  in  a  directed  network 
while  the  interdictor  tries  to  break  arcs  in  the  network  to  eliminate  all  possible  paths  for  the 
enemy.  Each  arc  (i,j)  has  associated  with  it  a  resource  expenditure  rij  which  is  required  to 
break  the  arc  and  the  interdictor  wishes  to  use  minimum  total  effort  to  disrupt  all  s-t  paths 
enemy  paths.  The  solution  to  the  problem  is  to  let  r/j  correspond  to  the  capacity  of  arc 
find  the  maximum  flow  from  s  to  <  subject  to  the  arc  capacities  and,  using  that  solution,  identify 
the  minimum  capacity  cut.  The  arcs  in  the  minimum  capacity  cut  are  those  which  should  be 
broken  to  eliminate  all  enemy  paths  using  minimum  total  effort. 

One  of  the  simplest  variations  on  the  above  problem  is  to  allow  only  a  limited  amount  of 
resource  to  break  arcs  so  as  to  leave  as  few  as  possible  arc-disjoint  paths  remaining  for  the  enemy 
to  use.  This  paper  will  be  primarily  concerned  with  a  generalization  of  this  problem: 

Problem  1.  An  enemy  wishes  to  move  as  much  of  a  single  commodity  from  node  s  to  node  i 
in  a  directed  network.  Each  arc  (i,j)  has  a  capacity  of  units  of  commodity  and  requires  of 
the  interdictor  an  expenditure  of  r,j  units  of  resource  to  break  the  aurc.  Partially  breaking  an 
arc  is  not  allowed.  The  problem  for  the  interdictor  is  to  minimize  the  maximum  amount  of  flow 
the  enemy  can  push  through  the  network  along  unbroken  arcs  where  the  enemy  is  constrained  by 
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the  arc  capacities  and  the  interdictor  is  constrained  to  use  no  more  than  the  R  units  of  resource 
available.  It  is  assumed  that  all  data  are  positive  integers. 

Problem  1  is  very  limited  as  stated  but  one  of  the  strengths  of  our  mathematical  programming 
approach  will  be  that  generalizations  are  easy.  Some  generalizations  that  will  be  considered 
include  multiple  resources  necessary  to  break  an  arc,  multiple  alternatives  for  breaking  an  arc 
using  distinct  resources,  multiple  sources  and  sinks,  multiple  commodities  which  are  needed  in 
specified  proportions,  and  the  ability  of  the  interdictor  to  partially  break  an  arc  using  a  fraction 
of  the  resource  needed  to  completely  break  the  arc. 

The  basic  problem,  with  minor  variations,  has  been  studied  before  in  [3-5].  (Related  work 
includes  [6-10]).  These  works  cite  military  applications,  not  drug  interdiction.  The  approach 
used  in  [3]  requires  the  network  to  be  source-sink  planar  so  that  the  network’s  dual  can  be  used, 
further  requires  the  enumeration  of  many  cuts  and  is  not  generalizable.  It  might  be  argued  that 
a  planarity  assumption  is  not  too  restrictive  in  practice,  although,  it  would  tend  to  rule  out 
adding  air  transportation  arcs  into  the  model.  However,  a  source-sink  planar  network  requires 
not  only  planarity  but  also  requires  that  the  source  and  sink  lie  on  the  outer  face  (exterior)  of  the 
network.  This  is  a  strong  assumption,  especially  when  a  drug  lab,  i.e.,  a  sink  node  for  precursor 
chemicals,  is  likely  to  be  found  hidden  in  the  middle  of  a  network  of  rivers  and  roads.  Further¬ 
more,  the  enumeration  of  cuts  is  an  exponential  process  in  the  worst  case  and  a  methodology 
based  on  this  enumeration  is  unsatisfactory.  The  methodology  used  in  [5]  is  limited  in  that  it  too 
requires  the  network  to  be  source-sink  planar  but  a  dynamic  programming  approach  is  employed 
which  avoids  enumerating  cuts  and  shows  that  the  source-sink  planar  problem  can  be  solved  in 
pseudo-polynomial  time.  Generalizations  as  simple  as  having  more  than  one  distinct  resource 
necessary  to  interdict  an  arc  would  be  impossible  except  in  very  simple  networks,  however.  (The 
difficulty  here  is  analogous  to  the  difficulty  which  arises  when  trying  to  generalize  the  dynamic 
programming  solution  of  a  knapsack  problem  to  multiple  resource  constraints,  i.e.,  exponential 
growth  in  complexity.)  The  methodology  described  in  [4]  does  not  require  planarity  and  a  spe¬ 
cialized  branch-and-bound  algorithm  is  developed  for  the  problem.  However,  that  approach  uses 
a  very  crude  bounding  procedure  which  would  be  difficult  or  impossible  to  generalize  to  multi¬ 
ple  resources,  multiple  commodities,  multiple  sources  and  sinks,  etc.  The  shortcomings  of  the 
aforementioned  approaches  invite  a  fresh  look  at  the  basic  network  interdiction  problem. 

In  the  following  sections  we  first  define  a  few  terms  and  notation,  consider  the  inherent  com¬ 
plexity  of  Problem  1  and  then  propose  an  integer  programming  model  whose  solution  we  claim 
yields  an  answer  to  the  problem.  We  then  prove  the  validity  of  the  model.  Next,  we  describe  a  few 
extensions  and  variations  on  the  basic  model  and  then  give  a  computational  example  previously 
found  in  the  literature.  Finally  we  consider  a  few  open  questions  and  make  concluding  remarks. 

2.  DEFINITIONS  AND  NOTATION 

G  —  {N,A)  will  denote  a  directed  network  with  node  set  N  and  arc  set  A.  We  will  usually 
refer  to  an  arc  as  an  ordered  pair  {i,j)  where  i,j  £  N,  although  we  can  also  refer  to  it  by  its 
number  k.  It  is  assumed  that  G  contains  no  self  loops,  i.e.,  no  arcs  of  the  form  (f,  i).  If  A'  C  A, 
then  G  —  A'  indicates  G  with  edges  A!  deleted  and  if  A'  C  A,  then  G  —  N'  denotes  G  with 
all  nodes  in  A'  deleted  along  with  all  arcs  incident  into  or  from  nodes  in  A'.  It  will  be  useful 
to  distinguish  two  nodes  s  and  t  with  s  ^  i.  Maximizing  flow  from  s  to  f  will  be  the  same  as 
maximizing  the  flow  along  an  extra  “return  arc”  {t,  s)  added  to  A. 

An  s-i  cutset  is  a  partition  of  A  into  two  subsets  N,  and  Aj  such  that  s  £  Aj  and  <  £  At. 
With  respect  to  that  cut,  an  arc  is  a  “forward”  arc  if  it  is  directed  from  a  node  in  A,  to  a  node 
in  At  and  it  is  “backward”  if  it  is  directed  from  a  node  in  T  to  a  node  in  S.  If  each  arc  (i,  j)  has 
a  capacity  «,;•  then  the  capacity  of  the  cut  is  the  sum  of  the  capacities  of  forward  arcs  associated 
with  the  cut. 

G  —  (A,  A)  may  also  denote  an  undirected  network.  For  an  undirected  network  the  definition 
of  a  cut  is  analogous  to  that  in  a  directed  network  and  the  capacity  of  the  cut  is  the  sum  of  the 
capacities  of  the  arcs  with  one  endpoint  in  Aj  and  the  other  endpoint  in  Nf  In  some  cases  it 
will  be  necessary  to  refer  to  an  undirected  arc  {i,j)  £  A  and  its  representation  as  a  directed  arc 
(i,  j)  and  a  directed  arc  {j,i)  in  anti-parallel.  Let  A'  denote  the  set  of  directed  arcs  (i,j)  such 
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that  i  <  j  and  (i,j)  €  A,  and  let  A"  denote  the  set  of  directed  arcs  (j,i)  such  that  (i,j)  €  A' . 
Thus,  each  arc  in  A  is  represented  once  in  A’  oriented  in  one  direction  and  is  represented  once 
in  A"  oriented  in  the  opposite  direction. 

In  proving  some  complexity  results  we  will  also  want  to  consider  an  undirected  graph  H  — 
{V,  E)  where  V  is  the  set  of  vertices  and  E  the  set  of  edges,  i.e.,  unordered  pairs  of  vertices.  The 
degree  of  a  vertex  v,  denoted  deg(u),  is  the  number  of  edges  of  the  form  (u,u).  A  self  loop  is  an 
edge  of  the  form  {v,v).  H  is  a  tree  if  it  is  connected  and  has  no  cycles.  If  E'  C  E,  then  H  —  E' 
denotes  H  with  all  edges  in  E'  deleted. 

3.  COMPLEXITY 

In  this  section  we  show  that  the  decision  problems  associated  with  Problem  1  and  two  variants 
are  NP-complete.  We  will  first  show  a  simple  transformation  of  the  binary  knapsack  problem 
(decision)  to  Problem  1  (decision),  which  are  defined  as  follows: 

Binary  Knapsack  Problem  (decision).  Given:  A  set  of  items  K  with  each  item  k  £  K  having 
a  positive  integer  profit  uj,  and  a  positive  integer  weight  rj,,  and  two  positive  integers  U'  and  R' . 
Question:  Does  there  exist  a  subset  K'  C  K  such  that  YlksK'  “i  ^  ^keK'  '"t  ^  ^ 

does  there  exist  a  set  of  items  whose  total  profit  is  at  least  U'  and  whose  total  weight  is  no  more 
than  R11 

Problem  1  (decision).  Given:  A  directed  graph  G  =  (JV,A)  with  distinguished  nodes  s  and  i, 
positive  integer  capacities  Uk  for  each  arc  k  £  A  and  positive  integer  resource  r*  required  for 
deletion  of  any  arc  k  £  A  and  two  positive  integers  U  and  R. 

Question:  Does  there  exist  a  subset  of  arcs  A'  Q  A  such  that  YlkeA'  ’’k  ^  R  And  the  maximum  s-t 
flow  in  G  —  A'  is  no  more  than  U,  i.e.,  does  there  exist  a  subset  of  arcs  whose  deletion  consumes 
no  more  than  R  units  of  resource  and  which  leaves  behind  a  network  with  maximum  s-t  flow  not 
exceeding  U? 

Theorem  1.  Problem  1  (decision)  is  NP-complete. 

Proof.  Consider  a  knapsack  problem  as  defined  above,  which  is  well-known  to  be  NP-complete 
(e.g.,  [11]).  Now  create  a  directed  network  G  =  (iV,  A)  with  two  nodes  s  and  t  and  for  each  item 
I:  G  A  in  the  knapsack  problem  create  an  arc  k  £  A  directed  from  s  to  f  with  capacity  =  uj, 
and  resource  requirement  rj,  =  r[.  Furthermore,  define  R  =  R'  and  U  =  —  U'.  Now, 

suppose  there  exists  a  subset  K'  C  K  such  that  uj,  >  U'  and  ’•i  <  -R'  Let  A' 

correspond  to  K' .  Then  it  follows  that,  because  of  the  simple  topology  of  the  network,  that  the 
maximum  flow  in  G  —  A'  is  at  most  Uh  —  U'  —  U  and  trivially  Ylkc.A'  ^  R-  Conversely, 
suppose  there  exists  a  set  of  arcs  A'  in  G  such  that  the  maximum  s-t  flow  in  G  —  A'  is  no  more 
than  U  and  YlkeA'  ’’k  ^  R-  Then,  letting  K'  in  the  knapsack  problem  correspond  to  A',  it  follows 
that  YlkeK'  “i  —  X]*  ^k  —  U  =  U'  and,  trivially,  J2keA'  ^'k  —  R-  Together  with  the  fact  that 
Problem  1  (decision)  is  cleMly  in  NP,  this  implies  that  Problem  1  is  NP-complete.  | 

Note  that  Problem  1  specialized  to  planar  networks  or  undirected  networks  is  stiU  NP-complete 
since  the  network  created  in  the  proof  is  planar  and  could  equally  well  have  been  undirected.  The 
proof  does  not,  however,  show  that  the  problem  is  NP-complete  in  the  strong  sense.  That  is,  the 
proof  leaves  open  the  possibility  of  a  pseudo-polynomial  time  algorithm  for  solving  Problem  1.  To 
show  that,  in  fact.  Problem  1  is  strongly  NP-complete,  we  show  that  the  following  specialization 
of  Problem  1  is  strongly  NP-complete. 

Problem  2.  This  is  the  same  as  Problem  1  except  that  interdiction  of  an  arc  requires  exactly 
one  unit  of  resource. 

So  in  Problem  2  the  interdictor  just  has  a  cardinality  constraint  on  the  number  of  arcs  he  can 
break  rather  than  a  general  resource  constraint.  The  proof  used  for  Theorem  1  does  not  follow 
through  for  Problem  2  since  a  knapsack  problem  in  which  each  item  weighs  exactly  one  unit  is 
trivial  to  solve.  To  prove  Problem  2  is  strongly  NP-complete,  we  state  the  problem  as  a  decision 
problem: 

Problem  2  (decision).  Given:  Directed  graph  G  =  (A,  A)  with  distinguished  nodes  s  and  t,  pos¬ 
itive  integer  capacities  of  Ujj  for  each  arc  (i,  j!)  £  A,  a  positive  integer  R,  and  positive  integer  U. 
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Question:  Does  there  exist  a  set  of  arcs  A'  with  \A'\  <  R  such  that  the  maximum  s-t  flow  in 
G  —  A'  has  value  {7? 

We  will  show  a  transformation  from  the  NP-complete  problem  “clique”  [11]: 

Problem  Clique  (decision).  Given:  An  undirected  graph  H  —  {y,E)  and  positive  integer  K. 
Question:  Does  there  exist  a  subgraph  of  H  which  is  a  clique  (complete  graph)  on  K  vertices? 

In  proving  that  Problem  2  is  NP-complete,  we  will  define  a  capacitated  directed  network 
with  respect  to  graph  H  and  show  that  there  exists  a  set  of  R  arcs  A'  in  G^  such  that  G^  —  A' 
has  a  maximum  flow  of  U  —  K  units  if  and  only  if  H  contains  a  clique  on  K  vertices.  R  will  be 
jj&l  —  (^) .  To  do  this,  we  prove  several  lermnas  which  use  an  intermediate  construct  illustrated 
in  Figure  1.  We  are  given  an  undirected  graph  H  =  {V,E),  which,  without  loss  of  generality,  has 
no  parallel  edges  and  no  self  loops.  From  H  we  create  a  capacitated,  directed  graph  G^  —  (N,  A) 
as  follows:  For  each  edge  e  G  E  create  a  node  ie  in  a  node  set  JVj  and  for  each  vertex  v  &  V 
create  a  node  j„  in  a  node  set  JV2.  In  addition,  create  special  nodes  s  and  t.  Now,  for  each  edge 
e  E  E,  direct  an  arc  in  G^  from  s  to  with  capacity  2  and  call  this  set  of  arcs  Ai.  Next,  for 
each  edge  e  =  («,  v)  direct  an  arc  in  G  with  capacity  1  from  to  and  direct  another  arc  with 
capacity  1  from  ie  to  Let  this  be  the  set  of  arcs  A2.  Finally,  for  each  vertex  v  eV  direct 
an  arc  with  capacity  1  from  jy  to  t  and  call  this  arc  set  Ag.  This  completes  the  construction  of 
G“  =  (JV,  A)  =  ({s}  U  {<}  U  All  U  7V2,  U  A2  U  A3). 


G"  =  {N,A) 


Figure  1.  Illustration  for  proofs  of  Lemmas  1  and  2. 


Lemma  1.  The  maximum  flow  in  G^  constructed  from  H  as  above  is  equal  to  the  number  of 
vertices  v  E  V  with  deg(v)  >  0. 

Proof.  In  other  words,  let  Vo  be  the  subset  of  V  containing  all  vertices  having  degree  0;  then, 
we  wish  to  show  that  the  maximum  flow  in  G^  is  |K  —  Vb|-  Suppose  that  H  is  a  tree.  Select  an 
edge  e  =  (ii,ti)  such  that  deg(ti)=l.  (Two  such  edges  must  always  exist  in  a  tree.)  One  unit  of 
flow  in  G^  can  be  routed  from  s  to  ie  to  j„  to  t.  Now  delete  e  and  v  from  H  to  create  another 
tree,  find  another  edge  e  =  (u,  v)  with  deg(i;)=  1,  route  another  unit  of  flow  from  s  to  the  new  ie 
to  the  new  to  t  and  repeat  until  H  consists  of  a  single  edge  and  its  incident  vertices.  At  each 
step  of  the  process,  it  is  possible  to  route  another  unit  of  flow  through  G^  because  the  path 
using  arcs  (s,ie),  (ie,jv)  and  {jv,t)  has  not  been  used  before.  Then,  with  one  edge  e  =  («,«) 
remaining,  we  can  route  one  unit  of  flow  from  s  to  ie  to  to  t,  and  because  the  capacity  of  (s,  ie) 
is  2,  we  ceui  also  route  one  unit  of  flow  from  s  to  ie  to  ju  to  t.  Since  a  tree  has  |P|  —  1  edges  we 
have  found  a  flow  of  |P|  units  from  s  to  <  and  this  is  obviously  maximum. 

Now,  consider  a  general  graph  H  =  {V,E).  Clearly,  the  maximum  flow  in  G^  is  at  most  |P— Vb| 
since,  if  deg(t))—  0,  it  is  not  possible  to  route  any  flow  through  j„  to  t.  Now  let  Hk  =  (14, Fi*) 
denote  a  connected  component  of  H  such  that  |£'ii;|  >  0,  i.e.,  deg(n)>  0  for  all  v  E  14-  Since  Hk 
contains  a  tree  we  can  route  one  unit  of  flow  (and  obviously  only  one)  through  each  node  jy  to  t 
for  each  vertex  v  E  14.  Thus  the  maximum  flow  in  H  equals  |  Uj,  14|  —  \V  —  Vb|.  I 
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Lemma  2.  Let  be  constructed  from  H  as  above.  Then,  there  exists  a  set  of  arcs  A\  C  Ai 
with  |A'i|  =  \E\  —  (^)  such  that  the  maximum  how  from  s  tot  in  —  A'^  is  K  if  and  only  if  G 
contains  a  clique  of  size  K. 

Proof.  (<=)  If  G  contains  a  clique  of  size  K,  let  E'  be  the  set  of  edges  not  in  the  clique  and 
let  Aj  =  {(s,ie)  I  e  G  E'}.  Of  necessity,  lAj|  =  IjE"!  =  |£1|  —  (^).  Then,  by  construction,  the 
maximum  flow  in  G^  —  A\  is  the  same  as  the  maximum  flow  in  G^~^  .  (If  (s,  ie)  created  from 
e  =  («,u)  is  deleted  from  G^ ,  we  might  as  well  delete  ie  and  edges  (ie,ju)  and  (ie,jv),  since 
no  flow  from  s  to  t  can  be  routed  through  them.  If  we  make  all  such  deletions,  G^  —  A'  is 
transformed  into  .)  But  by  Lemma  1,  the  maximum  flow  in  G^~^  is  K,  since  H  —  E' 

has  K  vertices  with  degree  /f  —  1  >  0  and  \V\  —  K  vertices  with  degree  0. 

(^)  Next,  suppose  that  it  is  possible  to  find  a  set  of  edges  A'j  C  Ai  with  \Ai\  —  |£i|  —  (^) 
such  that  the  maximum  s-t  flow  in  G  —  A\  is  K.  Let  E'  correspond  to  A'l  and  let  iV(  =  {»«  | 
ie  ^  El,  e  6  .  Then,  by  construction  G^  —  A'l  has  the  same  maximum  flow  as  G^  —  N[  and 

G^  —  N'l  =  G^~^' ■  Since  G^  —  N[  has  a  maximum  flow  of  K,  it  follows  that  H  —  E'  has  K 
vertices  v  with  deg(n)>  0  and  IPI  —  K  vertices  v  with  deg(n)  =  0.  Therefore,  H  —  E'  has  a  total 
of  (^)  edges  (none  of  which  are  in  parallel)  which  are  incident  to  exactly  K  vertices.  The  only 
way  this  can  occur  is  if  Jf  —  E'  consists  of  a  complete  graph  on  K  vertices  and  \V\  —  K  isolated 
vertices.  Thus,  H  contains  a  clique  of  size  K.  | 

Theorem  2.  Problem  2  is  NP-complete. 

Proof.  Problem  2  (decision)  is  clearly  in  NP.  For  any  instance  of  Problem  Clique  with  graph 
H  =  {V,  E),  we  will  create  a  capacitated  directed  network  G^  such  that  H  contains  a  clique  on  K 
vertices  if  and  only  if  G^  contains  a  set  of  arcs  A'  with  \A'\  —  R  such  that  G^—A'  has  a  maximum 
s-t  flow  of  K.  To  avoid  the  tedious  details  of  trivial  cases,  we  assume  that  \E\  >  1  and  K  >  2. 
Create  G^  by  first  creating  G^  from  H  as  described  before.  Then,  replace  each  arc  {ie,jv)  with 
|£1|  parallel  arcs  each  with  capacity  lj\E\  and  call  this  arc  set  A^.  Do  the  srune  thing  for  arcs  of  the 
form  {jv,t)  and  call  that  new  arc  set  As.  Then  G^  =  (N,  A)  =  ({s}u{f}UAri  GN2,Ai  Uj42UA3). 
Also,  let  .R=|£;|-(^). 

(■<=)  Let  E'  be  the  edges  of  E  not  in  the  clique  and  let  A'  =  {{s,ig)  \  e  €  E'}.  Then, 
and,  because  of  the  similarity  in  structure  between  G^  and  G^ ,  it  follows  from  Lemma  2  that 
the  maximum  flow  in  G^  —  A'  is  K. 

(=>•)  We  are  given  a  set  of  R  arcs  A'  such  that  the  maximum  flow  in  G^  —  A'  is  if.  If  A'  C  Ai, 
it  follows  immediately  from  Lemma  2  that  H  contains  a  clique  of  size  K,  because  of  the  similarity 
in  structure  between  G^  and  G^ ,  Now  suppose  that  A'  =  Aj  U  A2  U  A3  where  A2  U  A3  ^  0.  We 
will  show  that  the  maximum  flow  in  G^  —  A'  must  be  strictly  greater  than  K,  which  implies  that 
A'  cannot  contain  any  arcs  from  A2  or  A3  and  the  proof  will  be  complete.  Since  IA2  U  A3I  >  1,  it 
must  be  that  Ai  —  A)  contains  at  least  (^)  + 1  arcs  and  using  the  same  reasoning  as  in  Lemma  2, 
G^  —  A'l  has  a  maximum  flow  of  at  least  /f  +  1.  But  G^  —  A'  can  be  created  from  G^  —  A'j 
by  deleting  no  more  than  |JS|  —  (^)  <  |£;|  —  1  additional  arcs  taken  from  A2  U  A3.  Each  such 
deletion  from  G^  —  A'  reduces  the  maximum  flow  by  at  most  1/\E\,  since  the  capacities  of  these 
arcs  are  1/1E|.  Thus,  the  maximum  flow  in  G^  —  A'  is  no  less  than  A'  +  1  —  (1A|  —  1)(1/1A|)  = 
A'+1/|A|.  I 

Corollary  1.  Problems  1  and  2  are  NP-complete  in  the  strong  sense. 

Proof.  Problem  2  is  NP-complete  in  the  strong  sense  since  the  proof  of  its  NP-completeness 
did  not  require  any  entities  whose  number  was  dependent  on  the  numerical  values  Uij .  Problem  1 
is  therefore  NP-complete  in  tbe  strong  sense  because  it  can  be  specialized  to  Problem  2.  I 

The  proof  of  Theorem  2  uses  a  possibly  non-planar  graph  so  that  proof  does  not  imply  that 
Problem  2  is  NP-complete  for  planar  networks.  Indeed,  this  problem  cam  be  solved  in  polynomial 
time  for  “s-t  planar  graphs”  [6].  Likewise,  the  proof  does  not  imply  that  Problem  1  is  NP- 
complete  in  the  strong  sense  for  planar  networks  which  is  fortunate  since  a  pseudo-polynomial 
time  algorithm  for  this  problem  is  described  in  [5].  Finally,  we  note  that  if  Problem  2  is  further 
specialized  so  that  Uy  =  1  for  all  arcs  (i,j),  the  resulting  problem  is  of  polynomial  complexity; 
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Find  a  minimum  cardinality  cut  in  polynomial  time  using  a  maximum  flow  algorithm  and  then 
break  R  arcs  in  that  cut. 

4.  AN  INTEGER  PROGRAMMING  SOLUTION 

In  this  section,  we  give  an  integer  programming  model  which  will  solve  Problem  1.  We  first 
state  the  model  and  then  prove  its  correctness  starting  with  a  formal,  “min-max"  formulation  of 
Problem  1. 

The  formulation  of  the  model  is: 

Model  ID 

min  uij  A; , 

(iJ)eA 

s.t.  a,  -  aj  +  jiij  +  7y  >  0, 

at  —  a,>  1, 

'’•'i  Tij  < 
oi  e  {0,1}, 

AjiTy  G  {0, 1}, 

“D”  is  used  to  denote  that  this  model  is  for  a  directed  network.  Model  ID  is  based  on  a  modified 
dual  of  the  max  flow  linear  programming  formulation.  Essentially,  an  s-t  cut  is  identified  with 
all  £>,■  =  1  for  i  on  the  t  side  of  the  cut  and  a,-  =  0  for  all  i  on  the  s  side  of  the  cut.  The  value  of 
7,7  is  1  if  (i,;)  is  a  forward  arc  across  the  cut  which  is  to  be  broken;  A;  is  1  if  (i,  j)  is  a  forward 
arc  across  the  cut  but  it  is  not  to  be  broken;  and  all  other  Aj  Tij  are  0.  Thus,  we  see  that 
a  cut  is  identified  and  arcs  are  broken  in  that  cut  so  as  to  leave  as  little  remaining  capacity  as 
possible. 

In  order  to  prove  that  the  solution  of  Model  ID  solves  Problem  1  it  is  possible  to  proceed  from 
the  fact  that  the  max  flow  problem  is  totally  unimodular,  and,  conseqently  its  duetl  is  totMly 
unimodular  (e.g.,  [12]).  The  approach  taken  here  seems  more  direct,  however.  We  need  the 
following  lemma. 

Lemma  3.  The  dual  of  the  max  Sow  problem  has  an  optimal  solution  in  which  all  variables  are  0 
or  1. 

Proof.  The  max  flow  problem  is 


V  (i,j)  e  A, 

(1) 

(2) 

V  t  e  A, 


max  Xti, 

s.t.  ^  Xgj  —  ^  XJf  ~  Xts  ~  0, 
j  i 

~  =0, 
i  3 

xtj  -  xjt + = 0, 


Vi€N-{s,t}, 


0<a:y<Uj;,  V(i,j)eA, 

Xu  >  0, 

where  xts  corresponds  to  the  return  arc  which  has  been  added  to  the  network  going  from  t  to  s. 
The  dual  of  the  max  flow  problem  is 

min  Uij  $ij , 

(i.j)eA 

s.t.  Oj  -  aj  +  9ij  >0,  V  (i,  j)  6  A,  (3) 

“«-«j>l.  (4) 

0ij>0,  y(i,j)eA. 
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Let  {Ns,Nt)  correspond  to  the  minimum  capacity  cut  in  G.  Let  a,-  =  1  for  all  Nt,  let  a,-  =  0 
for  all  i  e  ,  let  6ij  =  1  for  all  arcs  (i,j)  which  are  forward  arcs  in  the  cut  and  let  all  other 
dij  =  0.  Clearly,  equation  (4)  is  satisfied  and  equations  (3)  can  be  seen  to  be  satisfied  by  checking 
against  the  four  classes  of  arcs 

(a)  ieN,,jeN., 

(b)  ieN.,j&N„ 

(c)  i  e  Nt,  j  G  N,,  and 

(d)  i  e  Nt,  jeNi. 

Thus,  this  solution  is  feasible.  Furthermore,  it  is  optimal  because  the  value  of  the  objective 
function  equals  the  capacity  of  this  cut  which  is  an  upper  bound  on  the  maximum  fiow,  but  the 
maximum  flow  equals  the  capacity  of  the  cut  by  linear  programming  duadity.  | 

We  can  now  prove  the  correctness  of  Model  ID. 

Theorem  3.  The  solution  to  Model  ID  solves  Problem  1. 

Proof.  A  model  to  solve  Problem  1  can  be  formally  stated  as  the  following  min-max  flow-based 
model. 

Model  2D 


min  max  Xt,, 

7€r  » 

s.t. 

>  j 

i  i 

=0, 

;■  i 

Xij  —  Uij  (1  —  ytj )  <  0, 

>  0, 


Vie  A-{s,t}, 


V(i,j)€  A, 

V(i,i)€  Au{(t,s)}, 


(5) 


where  T  =  {74  j  jij  €  {0,  l}V(i,  j)  6  A,  n,  <  R}- 

Now,  for  fixed  jij  the  dual  of  the  inner  maximization  problem  can  be  taken  giving  the  equivalent 
model: 

Model  3D 


s.t.  Of  —  Oj  -I-  dij  >  0, 
“«-“»>  1, 
a.  €{0,1}, 
Si;  €{0,1}, 


V(i,j)eA, 


Vie  A, 

v(i,i)eA. 


The  Oi  and  9ij  values  can  be  restricted  to  0  or  1  because,  for  fixed  values  of  yij,  the  inner 
maximization  of  Model  2D  is  Just  a  max  flow  problem  and  Lemma  1  applies.  Next,  Model  3D 
can  be  linearized  by  replacing  (1  —  7»j)  0i;  with  /?y  where  /Jy  e  {0, 1}  and  >  9ij  —  7y .  This 
yields: 

Model  4D 


min  ^  Uijfiij, 

s.t.  «(  —  Oj  +  9ij  >  0, 
a,  -  a,  >  1, 


v(i,i)eA, 
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Pij  T*i  ^ 

V(i,;)eA 

(6) 

(7) 

a,  G  {0,1}, 

WieN, 

e{o,i}. 

V  (i,i)  G  A, 

(8) 

Aj.T.y  e  {0,1}, 

V(i,i)GA. 

Next  we  argue  that  constraints  (6)  can  be  replaced  by  equalities.  Consider  any  optimal  solution 
in  which  Sy  =  0.  For  feasibility  and  7,j  can  be  either  0  or  1.  But  setting  to  1  would 
unnecessarily  increase  the  value  of  the  objective  function  so  it  must  be  0.  It  is  possible  that 
fij  =  1  in  some  optimal  solution  but  then  setting  y.j  to  0  creates  another  optimal  solution 
because  it  does  not  change  the  value  of  the  objective  function  and  maintains  feasibility.  Thus,  if 
6ij  =  0  it  may  be  assumed  that  Pij  —  jij  —  0,  or  equivalently,  ^ij  +  ytj  —  Oij  —  0.  Now  suppose 
that  in  an  optimal  solution  Oij  =  1.  In  that  case,  to  have  feasibility,  either  /Jy  =  1  or  yy  =  1  or 
both.  Suppose  that  Pij  ~  Then,  the  constraint  /?y  +  yy  —  5y  >  0  is  satisfied  and,  if  yy  =  1, 
there  exists  another  optimal  solution  in  which  yy  =  0.  Thus,  it  may  be  assumed  that,  when 
0y  =  1,  either  /3y  =  1  or  yy  =  1  but  not  both,  or  equivalently  /9y  +yy  —  0y  =  0.  Thus,  whether 
0y  =  1  or  Oij  =  0,  we  have  /?y  +  yy  —  0y  —  0,  which  implies  that  inequalities  of  constraints  (6) 
may  be  replaced  by  equalities. 

Since  /?y  +  yy  —  6y  =  0,  or  5y  =  /3y  +  yy ,  /3y  +  yy  can  be  substituted  everywhere  for  Sy . 
This  eliminates  constraints  (6)  and  (8)  and  yields  Model  ID.  | 

As  a  final  point,  note  that  constraint  (1)  can  be  eliminated  by  replacing  oj  by  1  and  a,  by  0 
wherever  they  appear. 


5.  EXTENSIONS  AND  VARIANTS 

In  this  section,  we  discuss  number  of  modifications  and  extensions  to  the  basic  model  which 
underscore  the  flexibility  of  the  mathematical  programming  approach. 

5.1.  Cardinality  Constraint 

Problem  2,  in  which  there  is  a  cardinality  constraint  on  the  number  of  arcs  interdicted  rather 
than  a  general  resource  constraint,  can  be  solved  by  solving  Model  ID  in  which  constraint  (2)  is 
replaced  by  yy  <  R.  We  will  refer  to  this  model  as  Model  IDC. 

5.S.  Partial  Arc  Interdiction 

In  this  scenario,  we  assume  that  by  applying  /yry  units  of  resource  to  arc  (i,  j)  where  0  <  /y 
<  1,  we  can  reduce  the  capacity  of  the  arc  to  (1  —  y .  Create  Model  2P  for  this  problem  which 
is  the  same  as  Model  2D  except  that  P  =  {yy  |  0  <  yy  <  lV(!,i)  G  A,  J2(i,j)sA  ^  ^}- 
Then,  making  transformations  analogous  to  those  used  to  develop  Model  ID,  we  obtain: 

Model  IP 

min  ^  Uijfiij, 

(>,3)eA 

s.t.  cti  -  aj  +  pij  +  yy  >  0, 

Oi  -  a,  >  1, 

'"•j  Ttj  < 

at  6  {0, 1}, 

0  <  yy  <  1 

0  <  Pij  <  1, 


V  (i,  j)  G  A, 


'dieN, 

V  {i,j)  G  A, 

V  {i,j)eA. 
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5.3.  Multiple  Sources  and  Sinks 

The  standard  way  to  handle  multiple  sources  and  sinks  would  be  to  create  a  super-source  s  and 
connect  it  to  the  individual  sources  s'  6  AT®  with  infinite  capacity,  unbreakable  arcs  and  to  create 
a  super-sink  t  and  connect  the  individual  sinks  i'  £  ,  to  it  with  infinite  capacity,  unbreakable 

arcs.  However,  the  fact  that  these  arcs  are  unbreakable  and  have  infinite  capacity  implies  that  in 
the  optimal  solution  the  super-source  will  be  on  the  same  side  of  the  optimal  cut  as  the  individual 
sources  and  the  super-sink  will  be  on  the  same  side  of  the  cut  as  the  individual  sinks.  Thus,  we 
know  that  a,i  =  0  for  all  s'  €  and  ati  =  1  for  all  t'  6  N'^ .  Therefore,  Model  ID  can  be 
modified  to  handle  multiple  sources  and  sinks  by  eliminating  constraint  (1)  and  substituting  a  0 
for  every  a,<  such  that  s'  £  IV®  and  substituting  a  1  for  every  af  such  that  t'  £  N^. 


5.4.  Undirected  Networks 

Next,  we  consider  the  analog  of  Problem  1  for  undirected  networks.  We  use  the  following 
proposition. 

Proposition  1 .  The  max-ffow  min-cut  theorem  holds  for  undirected  networks  as  well  as  directed 
networks. 

The  min  capacity  cut  identification  model  for  an  undirected  network  can  be  most  succinctly 
stated  using  A'  and  A"  as  defined  with  respect  to  an  undirected  network  G  =  {N,A)  in  Section  2; 


min  ^  UijOij, 

s.t.  Oi  —  Oj  +  6ij  >  0, 
s.t.  aj  —  ai  +  6ij  >  0, 
a,  -  a,  >  1. 


V(i,j)€A', 

V(j,j)6A', 


It  is  easy  to  verify  that,  associated  with  the  min  capacity  cut  (N,,Nt),  there  is  a  feasible  solution 
to  the  above  model  in  which  Oj  =  0  for  all  i  E  N,,  a,-  =  1  for  all  i  E  Nt,  6ij  =  1  if  i  £  Al,  and 
j  E  Nt  or  j  E  Ns  and  i  £  Nt,  and  otherwise  9ij  =  0.  Furthermore,  the  value  of  the  objective 
function  equals  the  capacity  of  the  cut.  That  this  solution  is  optimal  follows  from  Proposition  1, 
the  fact  that  the  dual  of  the  above  model  is  a  max  flow  model  for  an  undirected  network: 


maxxj,, 

X 

s.t.  ^  X,j  -  ^  Xjs  -  x„  =  0, 

j:(ii,j)eA'OA"  j:(j,s)£A'UA" 

=  0. 

3:(i,3)eA'uA"  3:(3,i)eA'uA" 

X]  H  =  0. 

3:(t,3)eA'uA"  3:(3,t)eA‘UA" 

Xij  +  Xji  <  Uij  , 


Vt  €  Ar-{s,<}, 


V(i,j)€A', 

V(i,j)£A'UA"u{(<,s)}, 


and  duality.  This  all  leads  to  the  analog  of  Lemma  1: 


Lemma  4.  The  dual  of  the  max  Sow  problem  for  undirected  networks  has  a  solution  in  which 
all  variables  are  0  or  1. 
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The  solution  to  Problem  1  for  undirected  networks  can  be  found  by  solving  the  following  min- 
max  model: 

Model  2U 


mini 

7€r 

max  xts , 

X 

s.t. 

E 

- 

E 

= 

0, 

j:(tJ)eA'uA" 

i.(j,s)eA'uA‘ 

' 

E 

Xij  - 

E 

= 

0. 

V 

i€N- 

-{s,0. 

j:(iJ)eA‘UA" 

j-(jA)^A'UA" 

E 

X,j  - 

E 

Xjt 

= 

0, 

j:(tJ)eA'uA" 

y.a.t)iA'uA‘ 

It 

Xi 

■j  +  X, 

a  -  “.7(1  - 

Hi) 

< 

0. 

V 

(».j)e 

A', 

Xij 

> 

0, 

V 

('■-i)  e 

yl'uA"u{(t,s)}, 

where  P  =  {yy  |  7y  6  {0,1}  V  (i,j)  6  A',  X)(ij)e4' ’’o ^  Following  a  similar  set  of 
transformations  and  using  Lemma  2,  Model  2U  is  seen  to  be  equivalent  to  the  simple  minimization 
problem: 

Model  lU 

min  ^  uyAj, 

s.t.  a{  -  aj  +  0ij  +  lij  >  0, 

«;■  -  Qii  +  I3ij  +  7ij  >  0, 

-  a.  >  1, 

Aj.Ty  e{0,i}, 


(9) 


V  i  6 

V  6  A‘. 


5.5.  Multiple  Resources 

Problems  with  multiple  resources  might  occur  in  two  cases.  In  the  one  case,  we  require  multiple 
resources  to  interdict  an  arc  where  interdiction  of  arc  {i,j)  requires  ryj  units  of  resource  I  for 
each  resource  I  €  L,  and  there  are  a  total  of  Rt  units  of  resource  I  available.  In  this  case,  in 
Model  ID,  equation  (2)  is  simply  replaced  by  the  set  of  equations 

X;  '••■i/T.y  V/6i.  (10) 

It  should  be  noted  that  if  some  of  the  ry;  values  were  negative  (suppose  by  interdicting  some  arcs, 
we  actually  capture  some  of  the  enemy’s  supplies)  the  substitution  used  in  going  from  Model  4D 
to  Model  ID  would  not  be  valid.  In  that  case,  it  would  be  necessary  to  use  a  model  of  the  form 
of  Model  4D  and  replace  equation  (7)  with  equations  (10). 

In  another  case,  it  might  be  that  the  interdictor  has  multiple  independent  resources  which  can 
be  used  to  interdict  an  arc.  For  instance,  three  different  types  of  aircraft  might  be  available  to 
attack  an  arc  and  exactly  one  of  the  types  of  aircraft  will  be  assigned  to  interdict  the  arc  or  none 
will.  For  this  case.  Model  ID  can  be  replaced  by 

min  ^  «yAj, 
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“1  +  Aj  +  7iii  >  0. 

; 

V(i,y)e 

at  -  a,  >  1, 

y~!  n}niji<Ri, 

'dlei, 

at  €  {0, 1}, 

V  t  6  iV, 

A;i  €  {0, 1}, 

V  (i,j)  e 

7iji  e  {0, 1}, 

V  (i,j)  e 

where  'j'yi  is  1  if  arc  (i,i)  is  interdicted  with  resource  /  and  is  0  otherwise,  and,  as  above,  Ri  is 
the  amount  of  resource  I  available. 


5.6.  MuUiple  Commodities 

Here  we  consider  a  situation  in  which  multiple  commodities  measured  in  the  same  units,  such 
as  kilograms,  are  being  sent  through  the  network  and  are  of  value  to  the  enemy  only  in  specified 
proportions.  For  instance,  suppose  it  requires  100  kilograms  of  chemicals  to  process  1000  kilo¬ 
grams  of  coca  leaves.  Then,  if  100  kilograms  of  chemicals  are  sent  through  the  network  along 
with  2000  kilograms  of  coca  leaves,  the  enemy  obtains  the  value  of  only  1000  kilograms  of  leaves. 

Let  bti  denote  the  number  of  units  of  commodity  k  which  the  enemy  requires  at  the  sink  node 
to  obtain  one  unit  of  benefit,  and  let  be  the  ^lmount  of  flow  of  commodity  k  across  arc  (i,j). 
The  min-max  problem  faced  by  the  interdictor  is  then 
Model  2MC 


minmajc  V, 

7€r  s,v 

s.t.  U  -  <  0, 

J  i 

i  > 

;■  i 

k 

^ijk  >  0, 


Vfe, 

VF, 

Vie  Ar-{s,t}, 


VF, 


V(i,i)6A, 

V  (f,i)  6  AU{Ut  {(**,«*)}}, 


where  T  =  {y,-,-  |  7,7  €  {0,l}V(i,ji)  G  A,  jij  <  /?}.  The  above  problem  definition 

allows  for  s  and  t  to  vary  by  commodity.  For  instance,  there  may  only  be  one  sink  node  t  where 
the  commodities  are  combined  but  the  sources  of  the  commodities  may  all  be  different.  If  all  the 
sources  were  the  same  and  all  the  sinks  were  the  same  the  problem  simplifies:  The  interdictor 
solves  the  single  commodity  problem,  the  enemy  finds  the  maximum  single  commodity  flow  in  the 
interdicted  network,  and  then  splits  that  flow  among  the  commodities  in  amounts  proportional 
to  ^k  for  each  commodity  fc'. 

Taking  the  dual  of  the  inner  maximization  problem  of  Model  2MC  yields 
Model  SMC 


7tl  a,  tf 

S.t.  Oik  -  ajk  -f  dij  >  0, 


V  k,  (ij)  €  A, 
VF, 
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IT*'*  ^ 

k 


>  0, 

V  Jb,  i  e  N, 

>  0. 

v(f,i)eA. 

Given  an  upper 

bound  9ij  on  9ij  allows  the  above  model  to  be  linearized  ; 

integer  problem; 
Model  4MC 

min  Uij0ij, 

s.t.  a.t  -  ajt  9ij  >  0, 

Vfc,  («,i)  6  A, 

Qftfcfc  -  >  0, 

VF, 

i^ij  +  7*j  “ 

v(«.i)eA, 

<  1. 

k 

^  T>j  ^ 

0J)eA 

Oij,  /?(;•  >0- 

V  (»,;•)€  A, 

Tij  6  {0,1}, 

V(f,i)€A, 

(Note  that  the  interpretations  of  the  a,  0  and  7  variables  given  for  Model  ID  are  not  valid  here 
since  these  variables  are  continuous.)  In  order  to  determine  a  suitable  value  for  6ij  for  a  given  7^ , 


note  that  is  the  dual  variable  on  the  constraint 


k 


in  Model  SMC.  If  the  right-hand  side  of  this  equation  were  increased  by  one  unit,  the  maximum 
change  in  the  objective  function  would  occur  if  this  allowed  one  additional  unit  of  commodity  k' 
to  flow  through  the  network  where  k'  =  argmax{6jC*}  and  this  allowed  t;  to  increase  by  units 
of  benefit.  Thus,  Model  4MC  is  valid  if  6ij  =  maxt{6J*}.  We  would  use  this  model  to  solve  the 
multicommodity  variant  of  Problem  1  since  there  is  no  analog  of  Model  ID  for  this  variant. 

6.  STRONGER  FORMULATIONS 

In  this  section,  we  describe  valid  inequalities  to  help  tighten  the  LP  relaxation  of  the  basic 
Models  ID  and  IDC  and  a  stronger  reformulation  of  Model  ID  based  on  cutsets.  (We  note  that  the 
reformulation  can  be  interpreted  as  introducing  additional  types  of  valid  inequalities  to  Model  ID 
but  it  suffices  to  view  the  process  as  a  reformulation.)  Initial  computational  experience  on 
Model  ID  [1]  indicates  that  fairly  large  problems  can  be  solved  without  the  need  to  introduce  cuts 
or  valid  inequalities.  However,  experience  has  shown  on  a  variety  of  combinatorial  optimization 
problems,  e.g.,  [13],  that  as  the  problems  become  larger  and  harder  to  solve,  valid  inequalities 
become  quite  useful.  Thus,  in  anticipation  of  solving  larger  problems,  we  discuss  two  types 
of  valid  inequalities.  The  first  type  of  inequality,  useful  for  Model  IDC,  is  derived  as  a  Chvatal 
inequality  [15]  with  a  strengthening  step  added.  The  second  type  of  inequality,  for  Model  ID,  is  a 
standard  knapsack  inequality,  e.g.,  [14].  The  usefulness  of  these  inequalities  and  the  reformulation 
is  demonstrated  in  the  next  section  which  gives  a  numerical  example. 


6.1.  Type  I  Valid  Inequalities  (Model  IDC) 

Consider  first  Model  IDC  which  has  a  cardinality  constraint  on  the  number  of  arcs  which 
can  be  broken  to  minimize  the  maximum  flow.  Consider  a  node  i'  and  a  maximum  set  of  arc 
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independent  paths  Pi'  from  i'  to  t.  Assume  that  \Pii\  >  R  and  let  A,'  be  the  set  of  arcs  on  the 
paths  Pii.  Then  the  following  “Type  I  inequality”  is  valid: 

Pij>\Pi'\-R-  (11) 


If  i'  is  on  the  t  side  of  the  optimal  cut  so  that  a;/  =  1  the  constraint  is  redundant.  However,  if 
i'  is  on  the  s  side  of  the  cut  so  that  a;-  =  0  the  constraint  simply  implies  that  at  least  |  arcs 
among  the  paths  will  cross  the  cut  and  since  at  most  R  of  those  crossing  arcs  can  be  broken, 
at  least  l/^/l  —  R  will  not  be  broken,  i.e.,  at  least  —  J?  of  the  jSij  for  {i,j)  on  those  paths 
must  have  value  1.  This  inequality  will  cut  away  a  fractional  solution  (i.e.,  it  is  violated  by  the 
solution  to  an  LP  relaxation  of  Model  IDC)  if  node  i'  can  be  chosen  so  that  op  <  1,  |P,/|  >  ii 
and  —  0  for  all  (i,  j)  G  A,-/ .  Such  a  situation,  if  it  exists,  can  be  found  in  polynomial  time  by 
finding  a  maximum  flow  from  i  to  f  in  G  where  the  capacities  of  arcs  are  0  if  /?.j  >  0  and  are 
1  if  =  0.  However,  there  may  also  exist  violated  inequalities  in  which  not  all  are  0  and 
these  would  be  harder  to  find;  an  integer  program  could  be  devised  to  find  them  but  this  might 
be  as  hard  as  the  original  problem.  It  is  theoretically  easy  to  generate  all  type  I  inequalities  a 
priori  using  a  maximum  flow  algorithm  coupled  with  an  enumeration  mechanism  which  would 
force  an  arc  to  be  in  a  maximum  flow  or  not.  However,  the  number  of  these  inequalities  may  be 
exponential  in  the  size  of  G  since  there  can  be  an  exponential  number  of  maximum  cardinality 
arc  independent  paths  from  any  node  i'  to  t. 

The  type  I  valid  inequality  can  be  derived  as  in  [15],  i.e.,  by  adding  multiples  of  constraints 
together  and  then  rounding  coefficients  to  integers,  although  we  add  a  strengthening  step  at  the 
end.  The  constraints  of  Model  IDC  can  be  written  as 


ai-aj  +  pij  +  jij 

a, 

at 

-  -in 

(‘J)eA 

at 

Pij  I  7ij 


>0, 

=  0, 

=  1, 

>  -R, 

V(bj)€A, 

(12) 

(13) 

€  {0, 1} 

V  i  6  IV, 

€  {0, 1} 

V(i,j)€A. 

As  above,  select  a  node  i'  and  a  set  of  paths  Pp  with  arc  set  A,i  and  such  that  |P,/|  >  R.  With 
multipliers  of  1  — c  for  e  sufficiently  small  and  positive  add  together  all  constraints  of  the  form  (12) 
for  (i,  j)  e  A,'  together  with  constraint  (13).  This  yields  the  redundant  constraint: 


(l-e)|/>w|«.-(l-e)|Pp|o,+  ^  (l-e)/?o-  ^  Hi  > -(I  -  ^)R- 


Noting  that  a/  =  1  then  yields 

(l-e)|P.-|a.v+  (1-OAi-  E  (l-0T<i  >(l-^)(l^i'|-«)- 

In  a  standard  fashion,  then,  all  coefficients  on  the  left-hand  side  can  be  rounded  up  to  integer 
values  while  maintriining  validity  of  the  inequality.  Then,  since  the  left-hand  side  is  integer  in  an 
optimal  solution  the  right-hand  side  can  be  rounded  up  to  the  nearest  integer  yielding: 

|P.-|a.-+  ^  Pii>\Pv\-R.  (14) 

If  Oil  =  0  in  an  optimal  solution  the  coefficient  on  Oji  is  irrelevant.  On  the  other  hand,  if  a;/  =  1 
constraint  (14)  is  redundant  and  remains  so  if  the  coefficient  on  a,'  is  decreased  to  IP,/  —  P|  which 
yields  the  type  I  inequality  (11).  This  inequality  is  clearly  stronger  than  (14). 


MCH  17:2-C 
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The  type  I  valid  inequality  can  be  extended  to  Model  ID  where  equations  (13)  are  replaced  by 

However,  in  forming  the  inequzJity,  a  good  choice  for  the  multiplier  on  constraint  (15)  is  not 
obvious.  Rather  than  pursuing  this  subject,  we  consider  completely  different  inequalities  for 
Model  ID. 

6.S.  Type  II  Valid  Inequalites  (Model  ID) 

Considering  constraint  (2)  alone  in  Model  ID,  it  is  clear  that  any  valid  inequality  for  a  knapsack 
problem  defined  on  that  constraint  is  valid  for  the  full  model.  Following  the  development  in  [12] 
suppose  that  yl'  C  v4  is  such  that  ry  >  R  but  r’y  <  R  for  any  A"  C  A' . 

Then,  .A'  is  a  minimal  independent  set.  The  extension  of  A' ,  denoted  E(A')  is  defined  by  E(A')  = 
A!  U  {(*',  j')  &  A  —  A'  \  Tiiji  >  Tij  V  {i,j)  6  A'}.  Then,  for  any  minimal  independent  set  A' ,  the 
following  inequality  is  valid  for  Model  ID: 

Yl  '••iTi;  ^  l^'l- 1-  (16) 

(•,;)€£U') 

We  refer  to  this  inequality  as  a  type  II  valid  inequality.  The  reader  is  referred  to  [12]  for  a 
discussion  of  other  types  of  knapsack  inequalities  which  could  be  applied  to  Model  ID. 

6.3.  A  Cutset-Based  Reformulation  (Model  ID) 

Here  we  discuss  a  cutset-based  reformulation  of  Model  ID.  If  all  cutsets  in  graph  G  were 
enumerated  and  there  were  not  too  many,  solution  of  Model  ID  would  not  be  trivial  but  it  would 
be  relatively  easy:  Let  A^  be  the  set  of  forward  arcs  in  cutset  C.  Solve  the  knapsack  problem 
on  each  cutset 


max  Y  “i;  Ti; . 

s.t.  Y  ’’•■f  I'i  ^ 

Tij  2  {O'  1}'  (I'i)  2  Ac, 

to  obtain  solution  7^  and  define  remaining  capacity  to  be  =  X)(i,j)6A'^(l  “  T«7)“ij-  Then, 
the  optimal  cutset  is  that  cutset  which  has  the  minimum  value  of  and  the  optimal  arcs  to 
break  are  those  arcs  (i,  j)  such  that  y-j  —  1. 

To  obtain  cutset-based  valid  inequalities  for  Model  ID,  however,  we  need  to  develop  a  cutset- 
based,  mathematical  programming  formulation  of  the  model.  Letting  C  denote  the  set  of  all  s-t 
cutsets  in  G,  the  following  formulation  is  probably  the  most  obvious: 

min  Y  “<jAj' 

(>.j)6A 

S.t.  xc  -  7ij  -  Pij  <  0, 

=  1, 
cec 

Y  '■■i  T.-j  <  R, 

(i,j)eA 

Tij -Aj  2  {0,1}, 

*c  6  {0,1}, 


V(i,j)2  A^, 


V  {i,j)  €  A, 

wcec, 


(17) 

(18) 
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where  fij  and  jiij  are  defined  as  before  and  Ep  =  1  if  C  is  the  optimal  cut  and  otherwise, 
xc  =  0.  The  difficulty  with  this  formulation  is  that  it  is  not  possible  to  write  down  a  simple, 
useful  relaxation  of  it.  If  constraint  (18)  is  deleted  the  optimal  solution  is  0.  If  con8traint(18)  is 
included  but  constraints  (17)  are  deleted  for  some  C,  the  optimal  solution  is  also  0.  The  following 
formulation,  called  Model  1C,  avoids  this  problem. 

Model  1C  uses  7,j  ,  and  a,-  defined  as  in  Model  ID.  With  respect  to  a  cutset  C,  let  Nf 
denote  the  nodes  on  the  s  side  of  the  cutset  and  let  denote  the  nodes  on  the  t  side  of  the 
cut.  New  binary  variables  if  and  if  are  defined  so  that  if  =  0  and  if  =  1  if  and  only  if  C 
is  the  optimal  cutset.  Finally,  let  ftp  denote  the  amount  of  resource  consumed  in  the  solution 
to  the  knapsack  problem  defined  with  respect  to  cutset  C,  let  be  the  set  of  arcs  {i,j)  such 
that  jij  =  1  in  that  knapsack  solution  and  let  be  the  set  of  arcs  (i,j)  such  that  =  1  in 
that  knapsack  solution.  The  reformulation  is 
Model  1C 

min  UijPij, 

s.t.  at  -aj  +yij  >  0, 

a.  <  0, 

<6  JVC 

;€JVC 

-  6?  -  T.i  <  0, 

<0, 

E  ’•'i  Tv  +  E(^  -  -  E(^  -  ^ 

(i,3)eA  cec  cec 

a,  =  0, 
at  =  1, 

All  variables  G  {0, 1}. 

Note  that  the  model  has  been  formulated  so  that  it  is  valid  even  if  C  is  replaced  by  some  subset 
of  C.  In  fact,  if  C  =  0  the  model  reverts  to  Model  ID. 

To  see  that  the  formulation  is  stronger  than  Model  ID  (we  do  not  prove  it  in  genereJ)  consider 
a  solution  (7*,/?*, a*)  to  the  LP  relaxation  of  Model  ID  in  which  aJl  variables  are  optimal  to 
the  integer  program  (IP)  except  one  and  the  paired  are  fractional,  i.e.,  the  LP  solution 
would  be  optimad  to  the  IP  if  7, were  forced  to  0  and  were  forced  to  1.  Such  solutions, 
in  which  the  optimal  cut  C  is  identified  but  not  proven  optimal,  can  and  do  occur  in  practice. 
In  one  of  these  LP  solutions  the  resource  constraint  will  be  tight  although  it  will  be  slack  in  the 
solution  to  the  IP.  Because  of  this,  if  Model  1C  is  formulated  with  C  =  {C*},  it  is  easy  to  see 
that  the  LP  solution  will  not  be  feasible  (the  modified  resource  constraint  will  be  violated)  for 
the  new  formulation  and  thus  Model  1C  is  stronger  than  Model  ID. 

7.  A  NUMERICAL  EXAMPLE 

In  this  section,  we  compaie  our  methodology  to  that  of  [4]  using  the  example  given  in  that 
paper.  The  test  network,  which  we  denote  G  =  (N,A),  is  given  in  Figure  2  and  the  capacity 
and  resource  necessary  to  break  each  arc  is  given  in  Table  1.  A  total  of  15  units  of  resource  is 
available  for  interdicting  arcs. 

G  is  undirected  and  it  has  multiple  sources  =  {1,2, 3, 4}  and  multiple  sinks  = 
{12,13,14}.  We  use  a  hybrid  of  Models  ID  and  lU  to  solve  this  problem  with  Oj  replaced 
by  0  for  eJl  i  6  A®  and  or,-  replaced  by  1  for  all  i  G 

Let  A  be  partitioned  into  Ast  and  Ast  where  Ast  denotes  those  arcs  incident  to  a  node  in 
AT®  or  and  Ast  denotes  the  complementary  set.  Arcs  in  Ast  niay  be  replaced  by  directed 


V(i,j)6A, 

VCGC, 

VGgC, 

v(i,i)GA:^,  cec, 

V(i,i)GA^,  CGC, 
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(j  indicates  3  is  a  source  node,  3  indicates  3  is  a  demand  node.) 
Figure  2.  Sample  network. 

Table  1.  Network  data. 


Arc 

Capacity 

Resource 

Arc 

Capacity 

Resource 

(1.5) 

60 

5 

(6,9) 

120 

4 

(1.8) 

70 

4 

(6,10) 

150 

6 

(1.6) 

60 

5 

(7,10) 

120 

6 

(2.8) 

50 

3 

(7.11) 

80 

4 

(2.6) 

50 

3 

(8,12) 

80 

4 

(2.7) 

60 

5 

(8,13) 

50 

5 

(3,6) 

100 

3 

(9,12) 

100 

5 

(3.7) 

80 

5 

(9,13) 

80 

4 

(4,6) 

50 

5 

(10,13) 

180 

6 

(4,7) 

100 

5 

(10,14) 

100 

4 

(4.11) 

80 

4 

(11,13) 

80 

5 

(5,8) 

60 

4 

(11,14) 

100 

6 

(5,9) 

60 

7 

arcs:  There  need  not  be  any  flow  into  a  source  node  so  any  arc  (i,  j)  with  i  €  may  be  replaced 
by  a  directed  arc  likewise,  any  arc  (i,j)  with  j  G  may  be  replaced  by  a  directed  arc 

(i,  j).  Let  A'gj,  be  defined  from  Ast  as  A'  is  defined  from  A. 

The  model  we  solve  then  is 


min  Uijpij, 

s.t.  a,-  -  aj  +  Pij  +  >  0, 

V  ((,  j)  €  A'^q<, 

Oi  -  aj  +  I3ij  +  >  0, 

V  (i,  j)  G  A'grp, 

otj  -ai  +  0ij  +  jij  >  0, 

V  {i,j)  G 

oii  =  0, 

V  i  G  N^, 

at  =  1, 

V  i  e  N'^, 

(i,j)eA^^VA's^ 

die  {0,1}, 

'iieN, 

A;,  7ii  e  {0,1)1/, 

V(i,i)  G^'. 

The  data  in  Table  1  yields  a  model  with  31  constraints  and  43  variables.  The  network  allows 
720  units  of  flow  if  no  arcs  are  interdicted.  The  optimal  solution  to  the  problem  was  obtained 
using  the  branch-and-bound  algorithm  within  LINDO  [16]  and  required  eight  branches  versus  the 
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sixteen  required  by  Ghare  ei  al.  The  initial  LP  relaxation  has  a  value  of  320  while  the  optimal 
solution,  obtained  on  the  third  branch,  has  340  units  of  flow  remaining.  The  interdicted  arcs  are 
(6,9),  (10, 13),  and  (10, 14),  while  the  unbroken  arcs  which  define  the  minimum  capacity  cut  are 
(1,8),  (1,5),  (2,5),  (4,11)  and  (7,11). 

As  asserted  in  the  previous  section,  the  solution  to  the  LP  relaxation  of  the  above  problem 
might  be  all  integer,  except  for  one  and  its  cohort  ~fij  being  fractional.  The  Oj  in  this  instance 
would  all  be  0  or  1  and  identify  a  cut  which  might  be  the  optimal  cut.  Such  is  the  case  for  this 
problem.  Adding  inequalities  of  type  II  and  reformulating  the  problem  based  on  the  identified 
cut  results  in  an  improvement  in  the  value  of  the  LP  relaxation  to  328.3.  The  branch-and-bound 
algorithm  then  requires  3  branches  to  solve  the  problem  optimally. 

To  test  inequalities  of  type  I,  the  resource  constraint  in  the  above  problem  was  replaced  by  a 
cardinality  constraint.  However,  the  solution  to  the  LP  relaxation  of  this  model  is  integral  and 
solves  the  IP  irrespective  of  what  integer  value  for  R  is  used.  Consequently,  it  was  necessary  to 
modify  the  until  a  fractional  LP  solution  weis  found;  a  value  of  if  =  3  was  used.  Evaluating 
the  LP  solution  yielded  a  type  I  inequality  based  on  four  arc-independent  paths  running  from 
node  9  to  the  sink  nodes.  The  LP  relaxation  of  the  model  proved  to  be  integer  optimal  after 
adding  this  one  inequality. 


8.  CONCLUSIONS 

This  paper  has  described  a  simple  network  interdiction  model  and  its  variants  in  which  an 
interdictor,  using  limited  resources,  interdicts  arcs  in  a  capacitated  network  so  as  to  minimize 
the  maximum  flow  that  can  be  pushed  through  the  network  by  an  adversary.  The  basic  problem 
is  shown  to  be  NP-complete  even  when  the  interdiction  of  an  arc  requires  exactly  one  unit  of 
resource.  A  new  integer  programming  model  is  developed  for  the  basic  problem  and  is  shown  to 
be  easily  modified  and  extended  to  handle  variants  and  generalizations  of  this  problem.  This  is 
in  contrast  to  the  methods  previously  known  which  cannot  handle  any  generalizations  such  as 
multiple  resource  constraints  or  are  restricted  to  planar  networks. 

In  anticipation  of  solving  large  network  interdiction  problems,  valid  inequalities  were  developed 
to  tighten  the  LP  relaxation  of  the  integer  programming  models.  A  cutset-based  reformulation  of 
the  problem  was  also  developed.  A  small  computational  example  illustrates  the  use  and  flexibility 
of  the  basic  model  and  shows  that  the  valid  inequalities  and  the  reformulation  can  achieve  quicker 
solutions. 
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